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Final  Progress  Report:  Mathematical  Imaging  and  Modeling  of 
Cortical  Spreading  Depression  and  Wound  Healing  Tom  Chou,  UCLA 


Foreword 

This  is  the  final  progress  report  for  “Mathematical  Imaging  and  Modeling  of  Cortical  Spreading 
Depression  and  Wound  Healing.”  The  research  conducted  under  this  grant  led  to  advances  in  image 
analysis,  modeling  cortical  spreading  depressions  waves  in  excitable  media,  and  inverse  problem 
theory.  It  has  resulted  in  ten  published  papers,  one  accepted  publication  and  two  other  submitted 
manuscripts. 

In  addition  to  the  originally  proposed  research,  new  but  related  directions  were  developed  in 
the  latter  stages  of  this  grant.  Specific  results  were  obtained  for  problems  involving  retinal  imaging, 
Bayesian  inference  of  stochastic  processes,  and  modeling  retinal  detachment.  These  results  built 
upon  results  and  techniques  developed  from  investigations  undertaken  in  the  earlier  part  of  the 
grant. 


List  of  Figures 

Figure  1:  Graph  cuts  algorithm  for  image  segmentation  and  shape  priors  for  object  identification. 

Figure  2:  Results  from  mathematical  models  of  corticol  spreading  depression  waves. 

Figure  3:  Schematic  of  delaminated  retinal  tissue  and  phase  diagram  for  the  stability  of  a  detached 
blister. 

Figure  4:  Simultaneous  reconstruction  of  bond  potentials  and  bond  mobility  from  trajectory 
data. 


Statement  of  Problem  Studied 

The  original  research  problem  was  to  to  develop  the  mathematical  tools  to  track  dynamics  of 
wavefronts  as  observed  in  cortical  spreading  depression  (CSD).  Accurate  tracking  of  the  depression 
front  was  then  be  used  to  inform  mathematical  modeling  of  the  biophysical  and  physiological 
mechanisms  that  govern  CSD.  Exploring  these  problems  also  led  to  applying  our  approach  to 
the  modeling  of  retinal  tissue  and  the  imaging  of  retinal  layers.  Machine  learning  approaches 
developed  for  front  tracking  were  also  extended  and  adapted  for  Bayesian  inference  of  simple 
stochastic  processes. 


Most  Important  Results 

Improved  segmentation  algorithms 

Our  work  on  imaging  and  tracking  moving  fronts  of  depressed  activity  across  cortical  tissue  has 
been  published:  J.  Chang,  K.  C.  Brennan,  and  T.  Chou,  “Tracking  monotonically  advancing 
boundaries  in  image  sequences  using  graph  cuts  and  recursive  kernel  shape  priors,”  IEEE  Trans¬ 
actions  on  Medical  Imaging ,  31,  1008  -  1020,  (2012).  In  this  approach,  we  exploited  ideas  Bayesian 
prediction  and  Kalman  filtering  to  compute  probability  distributions  for  propogating  fronts. 
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We  have  further  extended  this  method  by  reformulating  an  objective  function  for  statistical 
image  segmentation  so  that  faster  graph  cut  methods  can  be  applied.  We  applied  our  extended 
shape  prior  methods  common  images  as  shown  in  Fig.  1.  These  results  were  published:  Joshua  C. 
Chang  and  Tom  Chou,  “Iterative  graph  cuts  for  image  segmentation  with  a  nonlinear  statistical 
shape  prior,”  Journal  of  Mathematical  Imaging  and  Vision ,  49,  87-97,  (2014). 

Mathematical  modeling  vascular-activity  coupling  in  CSD 

A  manuscript  describing  the  effects  of  blood  vessel  dynamics  on  CSD  activity  has  been  submitted 
to  Biophysical  Journal.  In  the  model,  vascular  diameter  was  coupled  to  potassium  concentration. 
Changing  vascular  diameter  affects  the  flow  of  oxygenated  blood  to  the  tissue,  affecting  its  ATP 
production  rate,  and  ultimately  influencing  the  dynamics  of  the  recovery  of  the  depressed  activity 
tissue.  Figure  2  recapitulates  the  main  results  of  our  investigation: 

Modeling  Mitosis  and  Cellular  Cortical  Lamination 

This  study  provided  a  simple  method  for  estimation  of  cell-cycle  parameters  in  the  developing 
cerebellum  from  gross  histological  measurements  using  observed  global  constraints.  It  was  de¬ 
termined  that  symmetric  in-plane  mitoses  of  cells  in  the  cerebellum  are  the  most  important  for 
proper  cerebellar  foliation.  A  manuscript  by  Joshua  C.  Chang,  Mark  Leung,  Hamza  Gokozan, 
Patrick  Gygli,  Fay  Patsy  Catacutan,  Catherine  Czeisler,  and  Jos  J.  Otero  has  been  accepted  to 
The  Journal  of  Neuropathology  and  Experimental  Neurology. 

Modeling  glutamate  release  and  calcium  elevation  in  CSD 

In  this  combined  experimental  and  computational  study,  we  demonstrate  why  glutamate  release 
during  CSD  is  the  best  correlate  for  neuronal  compromise  by  its  coupling  with  calcium  dynamics 
inside  neurons.  This  study  resolves  the  conflict  between  Grafstein’s  K-mediated  theory  of  CSD 
propagation  and  Van-Harreveld’s  glutamate-related  theory.  We  characterize  three  different  wave- 
speed  regimes  based  on  the  amount  of  glutamate  released.  A  manuscript  on  this  work  has  been 
submitted. 

Mathematical  models  of  blister  formation 

We  first  studied  the  problem  of  retinal  detachment,  which  shares  many  mechanistic  similarities 
with  epithelial  blister  formation.  Our  results  are  summarized  in  the  Fig  3.  Our  main  result  the 
phase  diagram  that  maps  out  when  a  retinal  detachment  is  unstable,  or  is  stable  and  will  likely 
heal.  We  have  begun  to  generalize  this  problem  in  two  directions:  (1)  Developing  models  for  skin 
blisters,  and  (2)  Developing  a  more  detailed  understanding  of  the  tissue  attachment  condition. 
A  skin  blister  shares  many  similarities  with  retinal  blisters,  but  is  simpler  because  there  are  no 
known  active  fluid  pumps  and  fluid  does  not  pass  through  an  unruptured  blister.  However,  the  sub¬ 
blister  drainage  mechanism  is  more  complex,  involving  the  spatial  and  geometric  properties  of  both 
capillaries  and  lymph  vessels.  Moreover,  osmotic  pressure  from  fluid  that  leaks  out  of  capillaries 
is  expected  to  influence  the  ability  of  lymph  vessels  to  drain.  Finally,  we  wish  to  develop  a  more 
quantitative  understanding  of  the  forces  attaching  dead  epithelial  at  the  blister  edge.  This  contact 
line  of  a  blister  will  be  modeling  using  plate  theory  for  the  thin  tissue  layer,  including  the  effects 
of  tension  and  bending  rigidity  of  the  detached  epithelial  layer. 
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Bayesian  and  regularized  approaches  to  inverse  problems 

Finally,  new  methods  for  inverse  problems  were  developed.  This  line  of  research  arose  from  our 
work  in  quantitative  neuroimaging.  Here,  we  developed  path-integral  based  methods  for  recon¬ 
structing  coefficent  functions  in  PDEs  from  measured  data.  This  new  approach  incorporate  the 
data  mismatch  terms  into  an  “energy”  functional.  Treating  the  target  functions  to  be  reconstructed 
as  statistical  variable  themselves,  and  functionally  varying  the  data-infused  energy  functional  with 
respect  to  the  target  functions,  we  find  maximum  likelihood  solutions  the  target  functions.  How¬ 
ever,  the  energy  functional  also  includes  regularization  terms  reflecting  prior  knowledge  of  the 
expected  correlations  of  the  functions  to  be  reconstructed.  In  this  way,  we  can  reconstruct  prob¬ 
ability  distributions  over  the  target  functions,  providing  a  means  for  uncertainty  quantification. 
These  new  results  are  described  in  two  manunscripts,  one  is  published:  J.  Chang,  V.  Savage, 
and  T.  Chou,  “A  path  integral  approach  to  Bayesian  inference  for  inverse  problems  using  the 
semiclassical  approach,”  Journal  of  Statistical  Physics ,  157,  582-602,  (2014),  and  the  other  has 
been  submitted:  J.  Chang,  P.-W.  Fok,  and  T.  Chou,  “Simultaneous  Bayesian  reconstruction  of 
diffusivities  and  bond  potentials  using  path  integrals,”  Submitted  to:  Physical  Review  X,  (2014). 

In  this  second  project,  we  applied  our  functional  Bayesian  inverse  approach  to  directed  Brow¬ 
nian  motion  in  the  context  of  single-molecule  force  spectroscopy  (SMFS).  SMFS  is  a  ubiquitous 
experimental  method  used  to  determine  the  bond  potential  of  macromolecules  or  the  adhesion 
potential  of  macroscopic  materials.  By  repeatedly  applying  a  pulling  force  and  measuring  the 
bond  displacement,  data  in  the  form  of  “trajectories”  are  obtained.  These  trajectories  were  then 
used  to  simultaneously  reconstruct  both  the  bond  potential  and  the  bond  coordinate  mobility. 
This  approach  goes  beyond  the  current  state-of-the-art  that  uses  work  theorems  which  cannot  re¬ 
construct  bond  coordinate  mobilities.  Fig.  4  shows  the  simultaneous  reconstruction  of  a  complex, 
multi-minimum  bond  potential  and  the  bond  mobility.  Uncertainty  quantification  is  indicated  by 
the  yellow  bands.  Our  regularization  method  is  able  to  find  the  right  underlying  target  functions 
in  this  ill-posed  inverse  problem. 
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Figures 


Figure  1:  Left:  A  pixel-by-pixel  representation  of  the  graph  cuts  algorithm  for  deciding  which  pixels  belong 
to  foreground  or  background.  Segmentation  by  graph-cuts  is  computationally  more  efficient  than  other 
methods  such  as  finding  level  sets.  Right:  Implementing  the  graph  cuts  method  into  our  shape-prior 
algorithm,  we  are  able  to  discern  objects  in  the  presence  confounding  backgrounds. 
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Figure  2:  Left:  Shown  top  to  bottom  are  the  extracellular  potassium  level,  the  cellular  oxygen  concentration, 
the  potassium  ion  pump  flux,  and  the  typical  capillary  radius  for  a  typical  vascular  coupling.  Plots  are  time 
courses  at  a  fixed  position,  780/im  downstream  of  the  original  stimulus.  After  a  short  dilation  period,  the 
vessels  constrict  significantly  to  about  40%  of  their  rest  radii,  before  recovering.  Right:  Shown  (top  to 
bottom)  are  the  extracellular  potassium  concentration,  the  potassium  pump  flux,  and  the  cellular  oxygen 
level.  When  vasculature  is  coupled  to  extracellular  potassium  concentration,  the  extra  depletion  of  oxygen 
due  to  supply  constraints  causes  the  pump  to  operate  at  a  slower  rate  than  when  vasoconstriction  does  not 
occur.  As  a  consequence,  the  inward  potassium  flux  through  the  pump  is  decreased.  The  result  is  prolonged 
recovery  time.  The  vascular  coupled  model  (dashed)  is  contrasted  with  the  uncoupled  problem  (solid). 
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phase  diagram  for  blister  radius  Rx 


Figure  3:  Left:  Schematic  of  a  retinal  blister.  Ingredients  in  our  mechanical  model  include  vascular  pressure 
Pc,  subretinal  blister  pressure  ,  and  intraocular  pressure  Piop-  The  retinal  tension  T  is  balanced  by  the 
various  pressures  and  the  fluid  flow  (not  shown)  arising  from  the  retinal  pigment  epithelium  (RPE)  layer. 
Right:  Results  of  our  model  recapitulated  in  a  parameter  diagram  (the  effective  RPE  pump  flux  and  the 
effective  retinal  elasticity.  Domain  where  retinal  blister  are  stable,  or  unstable  to  growth  or  shrinkage,  are 
indicated. 


Figure  4:  Simultaneous  path-integral-baed  Bayesian  reconstruction  of  bond  potentials  and  mobilities.  The 
target  functions  are  indicated  by  the  red  dashed  curves.  As  more  trajectory  data  is  used,  the  most  likely 
reconstructions  approach  the  ground-truth  target  functions.  Moreover  the  95%  confidence  levels  are  found 
through  uncertainty  quantification. 
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